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Abstract. We use a new combinatorial technique to prove the optimal 
interior partial regularity result for LP -vector holds with integer fluxes that 
minimize the LP -energy. More precisely, we prove that the minimizing 
vector fields arc Holder continuous outside a set that is locally finite inside 
the domain. The results continue the program started in [PR], but this 
[ paper is self-contained. 



o 

Q 

03 



> 

On ■ Contents 
O . 

1. Introduction 2 

1.1. The result 2 

1.2. Using Calculus of Variations to construct U(l) -bundles 2 
psj ■ 1.3. Future steps: nonabelian bundles 4 

1.4. Relation to the regularity theory for harmonic maps and outline 

of the paper 5 

^ ' 2. The e-regularity theorem 7 

2.1. Smirnov's decomposition of 1 -dimensional normal currents 8 

2.2. Encoding the useful information in a graph 10 

2.3. Proof of the e-regularity 13 

2.4. A classical consequence: C°' a -regularity 19 

3. For minimizers X , weak convergence implies strong convergence 20 

3.1. Interpolant construction in the regular case 21 

3.2. Proof of Theorem 3.1 23 

4. The regularity result 25 

4.1. Dimension of the singular set 25 

4.2. Singular set of weak limits of minimzing vector fields 26 

4.3. Monotonicity and tangent maps 27 

4.4. Stationarity and dimension reduction for the singular set 28 

5. Stationarity and monotonicity 30 

5.1. Stationarity formula 30 

5.2. Monotonicity formula 31 
References 35 



O 

o 



2 



MIRCEA PETRACHE 



1. Introduction 

1.1. The result. In this work we consider vector fields X £ L P (B 3 , IR 3 ) such 
that for every a £ B 3 and for almost every r < dist(a, dB 3 ) , 

/ X-ueZ, (1.1) 

where v : dB 3 (a) — > S 2 is the outward unit normal vector. We call vector 
fields satisfying such flux conditions vector fields with integer fluxes. 

We observe that for p > 3/2 this class reduces to the divergence- free vector 
fields, and therefore we just look at the "interesting" range p £ [l,3/2[. It 
is clear that this class of vector fields is closed under LP -convergence. The 
compactness result for weak convergence holds only for p > 1 and is the main 
result of our recent work with Tristan Riviere [PR], and is based upon the 
introduction of a good distance between slices. 

Here we concentrate on the interior regularity of X £ L%(B 3 , IR 3 ) that are 
minima of the L p -energy ||X|| iP ( B 3). We say that X is a minimizer if it 
achieves the minimum in the following problem for some L p -function <fi defined 
on the boundary of B 3 and having integer degree, i.e. J QB3 (f> £ Z. 

inf | \X\ P : X £ L P Z (B 3 , M 3 ), u B s ■ X\ 9B z = J . (1.2) 

Note that without the constraint (1.1) the minimization in (1.2) yields the 
minimum X = regardless of the choice of <p. We recall the meaning of 
ubs ■ X\qb3 in Definition 1.4 and we show the existence of a minimizer for 
(1.2) in Lemma 1.6. This existence result depends on [PR] and the fact that 
in general it is nontrivial is implied by the properties of the trace present in 
[Pe3] . The main result of the present work is the following: 

Theorem 1.1. Let p e]l,3/2[, and let X £ L P (B 3 ,M. 3 ) be a minimizer. Then 
X is locally Holder-continuous away from a locally finite set EcB 3 . 

In the future work [Pe4] we will address the regularity up to the boundary. 

1.2. Using Calculus of Variations to construct 17(1) -bundles. In [PR] 
a first step was achieved towards the study of weak bundles with topologically 
nontrivial singularities. We recall here the main ideas of that approach. 
First of all, recall that U (l)-line bundles over 2-dimensional surfaces are classi- 
fied up to isomorphism by their first Chern class c\ , which for an 17(1) -bundle 
P over a compact surface E is expressible via Chern- Weil theory as 

ci(P)= f F A £ 2ttZ = /7 2 (S,Z), 

where Fa is the curvature of any connection A on P (see [Z]). By identifying 
the Lie algebra u(l) with M, we can identify Fa with an M-valued 2-form on 
S. In the "supercritical" dimension 3, a 2-form corresponds to a curvature 
if it assings integer volume to each closed surface (that integer corresponds to 
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the Ci of a line bundle restricted to the surface). By identifying 2-forms with 
1 -vectors in 3 dimensions, we arrive back at our definition of vector fields with 
integer fluxes. 

The idea started by [KR] was to study an energy minimization problem on 
bundles defined weakly as sketched above, the hope being that the minimizers 
would conserve some nontrivial information expressible in terms of C\ . 

Definition 1.2 ([KR]). Let Q C IR 3 be an open domain. We say that an 
LP -integrable 2 -form F on Q a curvature of a weak line bundle with group 
U(l) , if for all x £ Q and for almost all r > such that B(x,r) C Q, there 
holds 

[ i*F e z, 

JdB(x,r) 

where i : dB(x, r) — > M 3 is the inclusion map. We denote by jF%{£i) the class 
of such F . 

The above definition furnishes a definition of bundles in terms of slices on 
spheres, and a suggestive parallel can be made with the theory of scans as 
in [HR1],[HR2]. A first consequence of this parallel is the idea leading to the 
proof of weak compactness in [PR]. 

The main novelty of the above kind of definition in comparison with previous 
contributions, is that no assumption is made a priori regarding the existence of 
an underlying topological bundle structure, and one only assumes the existence 
of a curvature form respecting the Chern class constraint. This is the natural 
setting in which to study regularity problems, and in which to construct new 
bundles by minimizing the energy. In such way we "leave the minimizer free" 
to find the "most competitive" singularities, instead of forcing them on it at 
the beginning. 

Given an L p -integrable form <fi on dQ, consider the following minimization 
problem which is just the translation of (1.2) into the language of differential 
forms: 

inf{^|F|W/ 3 : FeJ%(n),? ga F = <l>y (1.3) 

The first positive result obtained in [PR] was the fact that J-f (Q) is closed 
under weak convergence. The precise statement of the result is as follows. 

Theorem 1.3 ([PR], Main Theorem). Suppose F n G ^f(^) converge weakly 
to a 2 -form F. Then F G 7f (fi) . 

The definition of the boundary value ig^F — is given in Section 5 of [Pe3]: 

Definition 1.4 ([Pe3]). Suppose F G J-"|(i? 3 ) , and suppose that <p is a 2- 
form in L p (dB 3 ) with j dB3 <f) G Z. Define the 2-forms F(p) := T*F , where 
T p : dB 3 B 3 is given by T p (a) = (1 - p)a for p G]0, 1[. We say that 
F G J-|^(-E> 3 ) if lim p _ s>0 + d(F(p), 0) = 0, where the distance d is as in [PR], 
[Pe3]. 
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It was shown in [Pe3] that 

Proposition 1.5. F n £ J% and F n — ^ F implies that ip n — > (ft for some 
LP -form (f) of integer degree, and that F £ , . 

We also recall that forms F correspond to vector fields X via the formula 
F P (U, V)=X p -(Ux V) for all U, V £ M 3 . 
The above results imply the existence of minimizers: 

Lemma 1.6. If <p is a 2 -form in LP [respectively, up to Hodge star duality 
with respect to the standard metric, an L p -function] on OB 3 having integer 
degree and the Definition 1.4 [respectively, its translation for vector fields] is 
used for the boundary value, then the minimum is achieved in Problem (1.3) 
[resp. (1.2)/. 

Proof. We give the proof in the language of forms. Consider a minimizing 
sequence Fj £ jF^AB 3 ) and extract a weakly convergent subsequence, which 
we label in the same way, abusing notation. We denote by F^ the limiting 
L p -form. From Theorem 1.3 we know that F^ £ F%(B 3 ) . Using Proposition 
1.5 we deduce that F^ £ J^JB 3 ). Thus Foo is the desired minimizer. □ 

Our Main Theorem 1.1 implies that a minimizer for the above problem, 
whose existence is proved above, actually gives an usual bundle (defined out- 
side isolated points) locally inside the domain. Therefore the program of con- 
structing weak bundles by variational methods works in the abelian case. 

1.3. Future steps: nonabelian bundles. The broader motivation for our 
work with abelian bundles is that they show some of the features of nonabelian 
(e.g. SU(2)-) bundles, without the complications due to the nonabelianity of 
the group. The next step after the proof of regularity present here, is indeed 
to attack the case SU(2) , that is the prototype of nonabelian bundles. In this 
case the classifying invariant is the second Chern class C2 £ LL 4 (M), defined 
for a 4-dimensional surface, again expressible in terms of curvatures, as 

c 2 (P) = / tr(F A AF A ) £ 8tt 2 Z. 

The above invariant is still present and significant for weak bundles, as shown 
in [U2]. The definition of weak bundles in a supercritical dimension (which in 
this case would be 5) by slicing [KR] can be stated also in this case (see also 
the treatment of [II], [12]): 

Definition 1.7. Let Qcl 5 be a domain. We call an L 2 -form F on Q with 
values in su{2) a representative of a weak SU(2) -bundle, if for all x £ Q and 
for almost all r > such that B(x, r) C f2 , there exists a gauge transformation 
g £ L 00 ^ 4 , SU(2)) and a W 1 ' 2 -connection A of an SU{2) -bundle over S 4 , 
such that the restriction Fb(t,x) of F to dB(x,r) ~ S 4 satisfies 

9~ l F B (r,x)g = dA + A A A. 
We call F%{Vt) the class of such F. 
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The idea is to look at the Yang-Mills energy j n \F\ 2 dH 5 in the above class, 
and the expectation is to obtain also in that case some analogy with the results 
for the £7(1) -case considered here. Definitions 1.2 and 1.7 coincide after replac- 
ing U(l),u(l),p by SU(2), su(2), 2, because in the abelian case g~ x Fg = F 
and A A A = , while in F = dA the regularity of such A directly follows 
from that of F after applying Fubini's theorem. 

Regularity results analogous to ours are not proved in the S77(2)-case, but 
a hint that they might be true comes from the singularity removal result of 
[TT]. Regarding the study of nonabelian bundles supercritical dimension, see 
also the more general works [Ti] and [DK], [DT]. 

The lack of an a priori given smooth structure will make the analogous of our 
result in the nonabelian case much different from previous works. The interest 
of the abelian case treated here is that it gives a new hint that the regularity 
result could be achievable. 

1.4. Relation to the regularity theory for harmonic maps and outline 
of the paper. Our regularity result parallels the following result of Schoen 
and Uhlenbeck (case p = 2), later extended by Hardt and Lin (for general 
p g]1,oo[) regarding minimizing harmonic maps. The result was proved for 
more general manifolds, but the special case stated here already presents the 
main difficulties. The more precise description the singularities is due to Brezis, 
Coron and Lieb. 

Theorem 1.8 ([SU],[HrL],[BCL]). Suppose u : B 3 -> S 2 is a map in W^B 3 , S 2 ) 
minimizing the L 2 -norm of its differential. Then u has Holder- continuous de- 
rivative outside a locally finite set EcB 3 . Moreover, u realizes a nontrivial 
degree around small spheres centered at each point in E . 

The analogy of our problem with the one of harmonic maps is also reflected 
by the fact that in our case the singularities also encode some topology, i.e. 
they all have a nontrivial degree. We decided however to prove such descrip- 
tion of the singularities in a future work, in order not to make the present 
article too heavy. 

The careful reader might be tempted to conjecture that our formulation of 
a minimization problem in terms of vector fields with integer fluxes can be 
reformulated in terms of harmonic maps. One could for example consider the 
minimization problem for the pullback u*usz of the volume 2-form on S 2 , 
via a map u G W 1,q (B 3 , S 2 ) . Since we are dealing with a 2-form u, the nat- 
ural regularity requirement for u corresponding to u*u G L p would then be 
u G W 1,2p . This is encouraged by the observation that the range of exponents 
corresponding to p g]1, 3/2 [ is q = 2p G [2, 3[ and gives precisely the Sobolev 
spaces for which the weak Jacobian d(u*u) of u is assured to be rectifiable 
and nontrivial. 

While at the level of function spaces there is no complication in sight, the 
problem is that the operation u i— > u*u is nonlinear in du, since u is a 2- 
form. Therefore there is no reason to think that the minimizing u should give 
a minimizing u*u, or vice versa. In general it is also not clear that -vector 
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fields are in bijective correspondence with forms u*u obtained as above from 
Sobolev functions u £ W 1,2p (B 3 , S 2 ) . In the (linear in du) case of 1 -forms, 
such representation is proved in [Pe2]. 

Our approach roughly follows the strategy of the regularity theory for har- 
monic maps. As in the harmonic map regularity proof, we derive and make 
use of a monotonicity formula and a stationarity formula (cfr [HrL] and [Pr] 
with our Section 5). In Section 2 we prove an e-regularity result, in Section 
3.1 we describe an analogous of the Luckhaus lemma [L], which helps showing 
the sequential compactness of minimizers. Then we proceed to the study of 
tangent maps and to the dimension reduction in Section 4. 

The techniques and results of sections Sections 2 and 3.1 are quite different 
from the approaches that we found in the literature, and might shed a different 
perspective also on the theory of harmonic map regularity. The main new 
observation is that the e-regularity can be studied on a simple model if we use 
the fact that the singularities come with an associated integer (the degree, or 
flux, of our vector field on small spheres surrounding the singularity). 
The structure that naturally arises is a weighted graph, having vertices that 
represent the singularities and edges representing the vector field's flow lines. 
Reducing to this model is allowed by the strong density result of Kessel, proved 
in [K] and summarized in [KR] . Its proof follows the strategy used by Bethuel 
to prove similar results for Sobolev maps into manifolds [Bl] [B2]. The precise 
result is the following: 

Theorem 1.9 ([K],[KR]). Suppose that Qcl 3 is an open set. Call TZ^{Q) 
the class of vector fields defined and smooth outside a finite set S = {a 1; . . . , a n } C 
Q, and having integer fluxes. Then 1Z^(VL) is dense in L%(Q) with respect to 
the L p - topology. 

The approximants to a minimizer (as given by Theorem 1.9) correspond 
then to normal 1 -dimensional currents. We are able to associate a weighted 
graph to vector fields in TZ^ , by applying a decomposition result of Smirnov 
[S] for normal 1 -dimensional currents (see Theorem 2.6). 

The e-regularity theorem is then obtained by a combinatorial reasoning on 
these graphs. It relies on an elementary minimax result (the famous "maxflow- 
mincut" theorem, [FF]). See the scheme (2.3) in the next section for a more 
precise overview of the proof. The same discretization method is the critical 
step also in the Luckhaus lemma, in Section 3. 

Acknowledgements I wish to thank Tristan Riviere for introducing me to 
the subject, for his encouragement, and for his many helpful comments and 
suggestions. I also thank Luca Martinazzi for some questions which helped 
clarifying the initial version of the paper. 
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2. The e -regularity theorem 



In this section our goal is to prove a so-called e-regularity theorem. This 
result states that if, for an energy minimizer X on a ball B the energy happens 
to be small enough, then X has no charges inside a smaller ball: 

Theorem 2.1 (e-regularity). There exists e p > such that for any minimizer 
X e L% of the L p -energy if Bf(x ) C B 3 and 



The main steps of the proof can be summarized as follows (see the scheme 



• We first approximate our vector field X E L% strongly in L p -norm by 
some smoother vector field X as in Theorem 1.9. 

• To X we associate a 1 -current T x in a classical way, and we apply 
to T x a decomposition result due to Smirnov [S] (see also the recent 
development [PS]). This result says that a normal current like T x can 
be decomposed via a measure (on Borel sets for the weak topology) 
fi x into a superposition of rectifiable integral currents supported on 
Lipschitz curves starting and ending on the boundary of T x . This 
result is described in Section 2.1. 

• Smirnov's decomposition fi x in our case (since the boundary dT x is 
supported on a discrete set) gives rise to a weighted directed graph 
G x , by grouping together the curves in the support of fi x with the 
same starting and ending point. These constructions are performed in 
Section 2.2. 

• We define a way of perturbing G x into another graph G' . For the 
underlying vector fields, this corresponds to perturbing X into a vector 
field X' that is (not smooth but) still in L%, and has energy bounded 
by the energy of X . We call these modifications elementary operations 
(see the definitions at the beginning of Section 2.2), and we use the same 
notation for operations on the graph G x and on the corresponding 
vector field X . 

• If X has little energy on a ball B , then we can perturb it by elementary 
operations into another vector field X' as above, and which has no 
charges inside B . This uses the classical "max flow/min cut" theorem 
on the graph G (see Section 2.3). 

• Finally, as the vector fields X approximate better and better the min- 
imizer X , since p > 1 we can apply the results of [PR] and extract a 
subsequence of the perturbed X' that converge weakly to a competitor 
for X . The comparison of X with the competitor gives a contradiction 
unless X has no charges in B, proving the result (see Section 2.3.2). 




(2.1) 



then 



divX = on 5 r / 2 (xo). 



(2.2) 



2.3 below). 
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directed graph 





(2.3) 



2.1. Smirnov's decomposition of 1 -dimensional normal currents. We 

build our constructions upon Smirnov's decomposition result for 1 -dimensional 
normal currents [S]. In order to state the results that we use, we need some 
preliminaries. 



Definition 2.2. A 1-current T in R 3 is called an elementary solenoid if there 
exists a 1 - Lip schitz function f : R — > R 3 with /(R) C spt(T), such that f,T 
satisfy 

T = V- lim -L/#[-T, A 
M(T) = I. 

In the spirit of the above definition, we can identify an oriented Lipschitz 
curve with a 1 -dimensional rectifiable current. We call Ci the set of all oriented 
curves of length < I , which we endow with the weak topology. All measures on 
paths described in this section will be positive, cx-finite measures, Borel with 
respect to the weak topology. The corresponding integrals are understood in 
the weak sense, i.e. 

S = Rdfi{R) is the current defined by 5(0) = / R(<j>)dn(R) for G V\R 2 

Definition 2.3. We say that a 1-current T is decomposed into currents lying 
in a set J C £> lj/oc (R 3 ) if there is a Borel measure \x supported on J such that 



T = J Rdfx(R), 
\T\\ = J\\R\\dn(R). 
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T G Ni ; ; oc (IR 3 ) is totally decomposed if the same \i also decomposes the bound- 
ary: 



dT = J dRdn(R), 

\\&r\\ = [ \\dR\\dfi(R). 



J .1 

Theorem 2.4. If I > , T G £>i(R 3 ),<9T = 0, then T can be decomposed 
into elements of Ci, with a measure \i of total mass M.(T)/£. Moreover, the 
following relations hold in the sense of measures: 

-\\T\\ > [ \\dR\\d^R), 



-||T|| = [ 5 b{R) dfi{R) = I 6 e ( R )dfj,(R), 
*■ Jc e Jc e 

where b(R),e(R) are the start and end points of R, respectively. 

We do not use the above theorem, but we cite it because in [S] it is the basic 
ingredient for the next two theorems, which we rely upon. Using Birkhoff's 
theorem (in the appropriate setting), Smirnov proves the following decompo- 
sition result. 

Theorem 2.5. T G X>i(lR 3 ), dT = 0, then T can be decomposed in elementary 
solenoids. 

For the case dT ^ there holds instead: 
Theorem 2.6. If T G Ni ) i oc (lR 3 ) then T can be decomposed as follows: 

t = p + g, 

imi = ii^ii + iiQii, 

dT = dQ,dP = 0. 

moreover Q can be totally decomposed into simple curves of finite length, i.e. 
into elements of := U^oC^. 

Remark 2.7. We now note some facts that follow easily from the constructions 
of Smirnov, but are not explicitly stated in his paper: 

(1) In the total decomposition of Q above, the paths have in general un- 
bounded (finite) lengths, but almost all of them (w.r.t. the decomposing 
measure fi) have b(R),e(R) on the support of dT = dQ . 

(2) If T corresponds to a regular vector field (i.e. for all test forms u , 
T(oj) = J u(X)dC 3 and X is regular), then the paths are composed of 
pieces of trajectories of the flow of X . 

(3) The functions b, e : — > M. 3 are continuous for the weak topology. In 
particular, given two Borel sets A,Bc M 3 , the set of paths 

{R : M(R) < oo,b(R) G A,e(R) G B} 

is Borel for the weak topology. 
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(4) Suppose that a 1 -current T decomposes via a measure \i on the space of 
1-currents. If a is a bounded Borel function on Vi(M. 3 ) , then v = afi 
induces by integration a 1 -current T a that is totally decomposed via 
\oc\fi, and satisfies 

T a = ±T and \\T a \\ < | |of| lioof^t) [ [T 1 ] |- 

Indeed, this is true for step functions a, and L 1 -convergence at the 
level of the decomposition induces weak convergence at the level of the 
decomposed currents. 

(5) The same result as above holds also in the case of a totally decomposed 
current T , with the analogous inequality holding also for the bound- 
aries: 

\\dT a \\ < \\a\\ L ^\\dT\\. 

2.2. Encoding the useful information in a graph. For vector fields X £ 
TZ^(Q) the decomposition of Smirnov allows to group the integral trajectories 
of X\_Q according to their start and end points: a generic trajectory could 
start or end on dQ or on one of the "charges" (i.e. singularities) of X . We 
encode this information in a weighted directed graph (i.e. a graph such that 
to each edge a positive number called "weight" and a direction are assigned). 
The weights in our encoding graphs keep track of how much of the flux of X is 
carried by each group of trajectories, and the direction of an edge encodes the 
direction of the corresponding trajectories. The grouping is done in such a way 
that there are no flux cancellations within the same group. Thus specifying 
the flux for a group of trajectories automatically gives a measure of the norm 
of the restriction of X to those trajectories. 

2.2.1. Elementary operations. The following kind of operations will be the ones 
that we perform on our encoding graphs: 

Definition 2.8. An elementary operation on a directed weighted graph G 
consists of multiplying by a factor a £ [—1,1] the weight of an edge, where 
multiplication of the weight by a negative factor a < means inverting the 
orientation and multiplying by \a\. 

We indicate by G -< G' the statement that G is achieved from G' , after 
applying finitely many elementary operations. 

We now define the elementary operations on the underlying X £ TZ^(Q). 
We use the same name because the two definitions correspond to each other 
in a natural way, as described in Section 2.2.3. 

Definition 2.9. Consider X £ TZ^(Q), which we identify with a current 
T = Tx as in Remark 2.7 (2), and to which we associate P, Q and a measure 
\i totally decomposing Q as in Theorem 2.6. An elementary operation on X 
consists in replacing X by the vector field corresponding to (Tx) a obtained as 
in Remark 2.7 (5), for some function a that only takes values in [—1, 1] and 
that is piecewise constant on a family of sets defined via b, e as in Remark 2. 7 
(3). 

We indicate by X <X' the property of X of being achievable after performing 
finitely many elementary operations starting from X' . 
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Remark 2.10. (1) It is immediate form Remark 2.7 (3) that X ^ X' 
implies \\X\\i,p < \\X'\\up with strict inequality unless \ot\ = 1 in all of 
our elementary operations. 

(2) is not invariant under elementary operations, since such opera- 
tions often create jumps in X . In general also the integer divergence 
condition is not preserved by these modifications. 

(3) From Remark 2.7 (5) it follows however, that for X £ D L P (Q) , 
any elementary operation sends X to a vector field X' £ L P (Q) having 
zero divergence away from the singular set of X . 



Figure 1. We represent schematically (i.e. we forget for a 
moment that we are in a 3 -dimensional setting, and we take 
Q to be a ball) the finitely many charges of our vector field 
X £ 7^00 n as black dots, and some of the supports of the 

rectifiable currents R of Definition 2.3, as thin lines. 

2.2.2. Grouping trajectories of X £ TZ^ n L 1 ^) . Consider X £ TZ^ n 
and the normal 1 -current Tx as in Remark 2.7 (2). 

Using Theorem 2.6, we can find a decomposition Tx = Px + Qx and a 
measure fix on := U^>cA that totally decomposes Qx into finite-length 
simple paths. 

Then note that, due to the special structure of X , d(Tx\-fl) is supported on 
dfl U {charges of X} . Also, by the total decomposition property of Qx , there 
holds 
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and b(R),e(R) G sptd(T x \-B) for fi x -a.e.R, so that we can decompose the 
set of finite length paths into disjoint Borel sets: 

n 

Coo = C U Cij, 
i,j=0 

where fi x (C) = and for all R G Cy there holds 

b(R) G A^,e(R) G Aj, 

where 

At := dQ n {sgn(X ■ v n ) = ±1} 

and 

-<4f ,i > enumerate the ± —charges of X, possibly with repetitions. 
By the decomposition theorem 2.6, if 



then 



j=0 i=0 

and for z > it is clear that fix(Cf) is equal to the charge of Af (see also 
Figure 2). 




: Cm 
'■ C02 
:Cu 

: C22 
'■ C10 
: C20 
: Coo 



Figure 2. In the example of Figure 1, we represent with differ- 
ent patterns the supports of curves belonging to different Cy 's. 
We omit the set if it has fix (Cij) = 0. 
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2.2.3. Associating a graph to a vector field. With the notations of the previous 
subsection, we associate to X the graph Gx (see Figure 3) which has the 
following features: 

• has vertices indexed by Af, i — 0, . . . , n, 

• has a directed edge A~ — > A+, for all < i, j < n, unless //(C y -) = 0, 

• any edge A~ — > A~j~ , it has weight nx{Cij) assigned to it. 




Figure 3. We superpose to the picture of Figure 2 the asso- 
ciated graph, where on top of each arrow we also describe its 
weight. The gray vertices Aq,Aq correspond respectively to 
start and end points of curves which lie on the boundary. 

Further, if G ^ Gx then we associate to G a vector field X <X such that 
G = G x , by the following procedure: 

• Fix a sequence Gx = Go ^ G\ -< ■ • • -< Gx = G such that G^+i is 
obtained from Gk by an elementary operation. We can still identify 
the vertices of Gk with those of Gx ■ 

• To each Gk we associate a function ak G L°°(fix), as follows. We start 
with cko = 1 ■ F° r k > if Gk+\ is obtained from Gk by multiplying the 
weight on A~ — > A~j~ by a G [—1,1] then we define ak+i ■= ctXCij a k + 
XCooXdjOtk- 

• Clearly ax G defines an elementary operation on X , and so 
we call X the vector field corresponding to (Tx)a N ■ 

2.3. Proof of the e-regularity. 

2.3.1. Modifications to eliminate charges in the regular case. In this subsection 
we restrict to vector fields X G TZ^(Q) satisfying the conditions of the e- 
regularity theorem, and we show that we can apply elementary operations 
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decreasing the energy while eliminating the charges of X. The main result is 
as follows. 

Proposition 2.11 (regular case). Suppose that X e IZ™ fl L p (Cl) and that 
Cl 3> Q is such that J dQ \X\ < 1 and J dn X • v = 0. T/ien i/iere exists a second 
vector field X e swc/i i/iai X ^ X and 

(1) X = X on tt\tt, 

(2) ll^llip(n) < 11^11^(0) an ^ 

(3) (divX)l_fi = 0. 

TTie inequality of 'point (2) is strict unless X already satisfies point (3). 

Proof. The main idea of the proof is to apply elementary operations to X , so 
that we cancel out the charges inside Q. Because of the above constructions, 
it is enough to do the corresponding operations on the graph Gx that encodes 
all the information that we need for the proof. 

Step 1: structure of the graph Gx- Consider the graph G := Gx defined 
in Section 2.2.3, and call 

• C + , C~ the sets of vertices of G corresponding to the interior charges 
of a given sign, 

• E 1 * 1 the sets of vertices of G corresponding to components of dQ with 
local charge ±, i.e. E 1 * 1 = {A^ }. 

The form of our graph is summarized in the following scheme, where we also 
indicate names for groups of arrows: 




The hypothesis j dn \X\ P < 1 implies that the arrows a ± have total weight 
less than 1 . This will be important in the sequel. 

Step 2: elimination of the singularities. We want to keep the arrows in 
cr ± fixed, and modify the other arrows via elementary operations so that the 
modified graph satisfies Kirchhoff 's law. This can be done as follows: 

• We keep (i.e. multiply by +1) all the edges which go directly from S + 
to S~ . Since these edges are not affected by the elementary operations 
done in the rest of the proof, we suppose from now on, without loss of 
generality, that there are no such edges. 

• Let's restrict to a connected component of our graph. Suppose first 
that it has the form drawn above (i.e. it is not degenerate): in this 
case we can find a maximal Kirchhoff subgraph K connecting S + to 
E~ , in the undirected graph 

£+ C- C + E~ . 

By the "max flow-min cut" theorem, after subtracting such directed 
subgraph, the remaining edges make a disconnected graph that has 4 
possible forms (where we keep the orientations as in the original G): 
(1) All arrows in v have been cut, but there are some edges joining E 
to some point charges. These charges correspond to singularities of 
X, for which at least 1/2 -charge flowed from/to E. In particular, 
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since the difference \a + \ — \ a~\ is constant during our construction, 
there must be an even number of such charges. This is not possible 
because the J" s \X ■ v^dV? was assumed to be smaller than 1. 

(2) The whole graph has been used, and we end up without leftover 
edges of the graph. Then again we see that J E \X ■ v^dH 2 is 
prohibited to be smaller than 1 , since in any charge connected to 
S 1 * 1 , the total wight of the arrows from/to the boundary dQ, is 
= | , and there are at least 2 such charges. 

(3) All arrows cr~ have been cut. Then also the arrows in a + have 
disappeared after eliminating the maximal flow, again because 
\a + \ — \a~\ is constant (equal to zero) during these modifications. 
Thus in this case all arrows outside v are canceled. Then we can 
multiply by zero the remaining arrows: these arrows are of posi- 
tive total weight since else we reduce to point (2), which is already 
excluded. Thus we strictly decrease the L p -norm of X . 

(4) The last case is the "generic" one: it could be that after the cut 
we are left with a graph of the form 

C~ C7+ 

£+ """" Tr 
C~ - C + 

It is shown in Lemma 2.13 that in this case it is possible to find 
another minimal cut that gives a graph as in (3) or in (2) instead, 
and we conclude the proof. 
The conclusion of this enumeration is that the only possible cases that 
allow any singularity at all inside Q and are compatible with the small 
boundary energy are the ones corresponding to the case (3) above. 
Observe that in this case we are sure to have canceled some edges i.e. 
we have decreased the energy of X , as wanted. 

□ 

Example 2.12. Consider a regular vector field in TZ,^ n L P (Q) that has 5 
singularities, one point having charge 1 a second point having charge 2, and 
the remaining points having charge —1 each (see Figure 4)- Suppose that the 
weights of the edges of the associated graph are as in Figure 4- We assume 
from the beginning that fix gives no weight to the curves that both start and 
end point on the boundary (such curves are anyways not affected by our ma- 
nipulations). The maximal flow showed on the right corresponds to any of the 
3 minimal cuts on the left. In general, no uniqueness of either the maximal 
flow or the minimal cut is guaranteed. In Figure 5 it is shown what happens 
next, in our manipulations. Once we fix the maximal flow of Figure 4, we 
change by elementary operations the flow lines of X , ending up with the graph 
on the left of Figure 5. Since this represents a flow, i.e. obeys Kirchhoff's law, 
the curves representing the modified vector field X are concatenated, i.e. that 
they all start and end on the boundary. This concatenation is u automatically 
done" by Smirnov's decomposition, since the associated current T x is totally 
decomposed (see Definition 2.3). The u canceled flow" on the right of the figure, 
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Figure 4. (A): a graph corresponding to a possible vector field 
X having 6 charges (of which the one represented by a larger 
circle is a double one). In the unoriented graph (B), we represent 
by dashed lines two minimal cuts separating the vertices with 
dashed boundaries; a non-minimal cut is represented by a dotted 
line. Observe that the flow through each of the 3 cuts in (A) 
is the same, but in (B) the sum of edge capacities is larger for 
the dotted line. In (C) we show the unique maximal obtained 
on (B) between the gray vertices. 



gives a measure of the amount of LP -norm of X gained this way. 

We must point out that the LP -energy improvement in passing from X to X 
depends also from factors not captured by the graph Gx itself, namely on the 
lengths and concentrations of the curves decomposing the associated current 
T x ■ But for our purposes a subtler analysis along these lines is not needed. 

Lemma 2.13. Under the hypotheses of Proposition 2.11 on X , suppose that 
a connected component of the associated graph Gx has the form 

x C- ,...« 

b"* (■ _ — c+ f 

where a minimal cut is given by the arrows in b, c,d,e. Then another minimal 
cut is given by the arrows in a,b. 



Proof. The fact that a, b give a cut is clear from the above diagram. We must 
show that such cut is a minimal one. 
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Figure 5. Continuing with the example of Figure 4, we show 
schematically on the left what remains after the cancellation of 
the charges (in terms of the associated graphs the three arrows 
of weight | actually are substituted by just one arrow of weight 
|, but we drew the picture to suggest that a procedure of "con- 
catenating arrows" is actually underlying the operation). On 
the right we show the flow that results after removing the max- 
imal flow graph out of the initial graph. In our charge removal 
procedure, we diminish the weights of our graph by the amounts 
in the right picture, so in this particular X has a smaller energy 
than X . 

We indicate by |x| the total flow through the arrows of the group labelled by 
x. First of all observe that by the zero total flux and small boundary energy 
hypotheses on X , 

M + l&l = |e| + |/|<i 

therefore, being b, c,d,e a minimal cut, by comparison with the above cut we 
obtain 

\b\ + \c\ + \d\ + \e\ < -. 

This implies that the total number of charges contributing to the vertices C + 
is the same as the number of charges contributing to C~ , and similarly for 
C + ,C~ . Indeed, suppose for contradiction that the numbers of charges con- 
tributing to C + ,C~ were not equal. Then the total flow \a\ + |c| + \d\ + |e| 
would be > 1, and this would contradict the fact that \a\ and |c| + \d\ + |e| 
are both < | . 

By the consideration in italics above, we obtain that 

\a\ + \d\ = \c\ + \e\, \b\ + \c\ = \d\ + |/|. 
Therefore, by definition of a minimal cut 

\b\ + |c| + |d| + |e| < \a\ + \b\ 
and this gives, using the previous computations, 

H > \c\ + \d\ + \e\ = \a\ + 2\d\, 
so \d\ = and the above inequalities are actually equalities, as wanted. □ 

2.3.2. The proof of e -regularity. 

Proof of Theorem 2.1. First of all, we may reduce to the case where the e- 
regular ball B(xo,r) of the theorem is the unit ball B = -8(0, 1), since the 
estimates and the function spaces considered are invariant under homotheties 
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and translations of M 3 . We call Cl the image of the initial B\ under this 
transformation. 

Step 1: fixing a small energy sphere. We claim that, for any small 

e > 0, we can find a positive measure set of radiuses p > 1/2 such that 
J dB \X\ P < 2e p . Indeed, if the opposite estimate would hold for a.e. p > 1/2, 
then we would obtain 



B 



\X\ P > I I \X\ p >e p , 



therefore 



/ \X\ p dH 3 > (1 -e)e , 
Jb 



and this contradicts our assumption for e small enough. Now from the above 
boundary energy bound by e p we get via Holder's inequality the following 
bound 

/ |X|<(2e p )^ 2 (S 2 )]^, 

J dB p 

and we choose e p such that the right hand side is equal to 1 . This gives 
the small boundary energy condition as in Proposition 2.11, and the zero 
flux condition follows from the definition of L%(B) and from the inequality 
\X-"b p \ < \X\. 

Step 2: passing to the approximants. We know that there exist Xk G 
TZ^(Cl) that converge to X in L p -norm. From the construction leading to this 
approximation it is clear that we can also further impose the convergence 

Xk\dB p — > X\g Bp , 

therefore for k large enough, Xk satisfies the properties required in Proposition 
2.11. Applying this proposition, we thus obtain Xk G L%(Cl) which are equal 
to X k outside B p and satisfy H-XjfcHx^ft) < H^fclLprn) (with strict inequality 
if (divX k )\-B p ^ 0) and (divX k )\-B p = 0. 

Step 3: a divergence-free competitor. By weak compactness of L^fl) it 
follows that a subsequence of the Xy converges weakly to some X G 
The zero divergence condition passes to weak limits, so divX = on B p . By 
sequential weak lowersemicontinuity of the norm, we also deduce 

||X|| p < liminf ||Xfc/|| p < liminf ||Xfc|| p = 

Since X was a minimizer, all the above inequalities must actually be equali- 
ties. We also observe that since the sequence Xk> converges both weakly and 
in norm, it must converge also strongly, to X . By examining the definition of 
elementary operations we also observe that the inequality |Xfc|(x) < 
holds almost everywhere for all k, and from it and the a.e. convergence it 
follows that the same inequality holds also in the limit. Since both X and X 
are minimizers it further follows that = |AT|(x) almost everywhere. 



Step 4: X is also divergence-free. We use the classical regularity theory, 
namely Lemma 2.15 (which applies since divX = 0) and Proposition 2.16 to 
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deduce that X is Holder-continuous in the interior of B r i 2 . It then follows 
that also divX = on B r / 2 , since in this case X £ L^(B r j 2 ). Indeed, 
using Theorem 1.9 it follows that X can be approximated by vector fields in 
TZ^(B r / 2 ) in the strong norm in L q for q > 3/2. But for such exponents 
the vector fields in n L q (Q) are all smooth (and in particular divergence- 
free, since the divergence is concentrated at their singular points). Thus by 
approximation also X is divergence-free. □ 

2.4. A classical consequence: C 0,a -regularity. From Theorem 2.1, using 
an extension by Peter Tolksdorff (and Christoph Hamburger) of the regularity 
theory first developed by Karen Uhlenbeck, it is relatively straightforward to 
prove the following extension of it: 

Theorem 2.14 (Holder version of the e -regularity). If X £ L| is a minimizer 
then we can find an e p > such that if on B 3 (xq) C B 3 the vector field 
X satisfies (2.1) then on B r / 2 (xo) the vector field X is a -Holder, with a 
depending only on p and with the Holder constant of X\b v/2 depending only 
on p and on | \X\ Ilp^) ■ 

In order to prove the above theorem, we use the conclusion that divX = 
of Theorem 2.1 and the Euler equation of the functional J n \X\ P to reduce to 
the by now classical regularity result for systems of equations due to the above 
cited authors. The main heuristic idea in play here is that roughly "divX = 
implies that X = V/ for some -function /". 

In order to use this idea while still keeping rigorous, we use the formulation of 
our minimization problem in terms of differential 2 -forms u instead of vector 
fields X. 

Lemma 2.15. The condition that a vector field X £ L%(Q) minimizes the LP - 
energy and satisfies divX = implies that the associated 2 -form uj £ ^f(^) 
satisfies locally in the sense of the distributions the following equations: 

j du = 

\ 5 (\u\ p - 2 u) = 0. 

Proof. The first equation is a trivial translation of divX = in our new setting. 
The second one is the Euler equation, and can be directly obtained from the 
requirement that u be minimizing, by using the perturbations u H- u + ed(p, 
for £ C^°(A 2 Q) and taking the derivative in e at e = 0. Since dip is exact, 
it easily follows that the perturbed form is still in J-|(0). □ 

With the result of the above lemma, we are exactly in the setting of [Ul], 
except that that article treats the case p > 2, while we are interested in the 
case 1 < p < 3/2. 

Luckily, the result of [Ul] was extended in [Tol], [To2] to the case 1 < p < 2. 
The article of Tolksdorf considers only the "basic case" where the equations 
concern a differential of a function instead of the generalization of exact differ- 
ential forms described by Uhlenbeck, but the setting in which Tolksdorf proves 
regularity can be translated without much effort into the one of Uhlenbeck, 
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and the techniques present there are not affected by the translation. Later 
on Hamburger [H] partially recovers the approach of Uhlenbeck for the whole 
range of exponents of Tolksdorf, but this work deals with homological mini- 
mizers (i.e. minimizers with respect to perturbations as in the proof of Lemma 
2.15 that keep the comparison forms in our class J-f(f2)) instead of just using 
the Euler equations as the other two works. In particular, the range of expo- 
nents p g]1, 2[ is recovered from the range p > 2 via a duality argument where 
the requirement of dealing with minimizers is involved. Due to these consid- 
erations, we can safely state the following version of these regularity results. 
A recent treatment of the regularity for differential forms including the case 
needed here is [BS]. 

Proposition 2.16 ([Tol],[H],[BS]). If u e L p (A 2 tt) satisfies the equations of 
Lemma 2.15 in the weak sense, then u is a -Holder, with a depending only on 
p and with the local Holder constant of X\b t/2 depending only on p and on 
\\X\\up(b t ) f or an U ball contained in Q. 

From the above lemma and the proposition, it is straightforward that The- 
orem 2.14 holds. 

3. For minimizers X, weak convergence implies strong 

CONVERGENCE 

In this section we prove the following compactness result: 

Theorem 3.1. Suppose Xf. G L%(B) are minimizers of the L p -energy, and 
that Xk — X weakly in L p . Then X is also a minimizer and Xk —t X 
also IP -strongly on any ball B(0,r),r < 1. In particular, any sequence of 
minimizers of bounded energy has a strongly convergent subsequence. 

It is a classical result that strong convergence can fail while weak conver- 
gence holds, only if some energy is lost in the limit. Thus, it remains to prove 
that the energy of X on B r is not lower than the limit of the energies of 
the Xk on the same ball. The fact that any X obtained as a strong limit of 
minimizers is a minimizer itself follows from the strong local convergence. 

The idea of the proof is to introduce a small parameter e > and to con- 
struct an interpolant Xk G L%(B) that equals Xk on B \ B r+e and X inside 
B r , in such a way that the energy of Xk in the small spherical shell B \ B r+e 
goes to zero as e — > . This allows us, using the minimization property of Xk , 
to bound from above the energy of Xk on B r , by the energy of X on the same 
ball. 

For the proof of the e-regularity, it is enough to be able to do the construc- 
tions for vector fields in TZoq. The interpolation construction faces again a 
problem related to possibility that (the approximant of ) X — Xk have some 
singularities in the small shell B r+e \ B r . We deal with this situation again by 
choosing shells where on the boundaries X — Xk does not have large energy 
for k large, and by applying the singularity removal operations of Proposition 
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2.11 from the e-regularity proof. After these elementary operations, we are 
reduced to an easier situation (see Figure 6), where the curves of the Smirnov 
decomposition of our vector fields all move from one boundary of the shell to 
the other. In this simpler case, the interpolation can be done via an auxiliary 
function / satisfying a Neumann boundary value problem in the shell, and 
the scaling of the classical energy bounds as the thickness of the shell vanishes 
(see Lemma ref scaling), are strong enough for our purposes. 

3.1. Interpolant construction in the regular case. The result on the 
existence of the interpolants that we need is the following. 

Proposition 3.2. There exists a constant C depending only on our exponent 
p from above, such that the following holds. For any numbers R and e such 
that R > 1 + e > 1 , for any Y £ IZoo^Br) having zero flux through dB\ +e 
and through dBi , having no singularities lying on these two boundaries, and 
satisfying 

[ \Y\dU 2 < i 

JdB r z 

for r = 1 and for r — 1 + e, there exists another vector field Y £ Hoo{Br) , 
such that 

• Y = Y on B X) 

• Y = outside B\ +e , 

— ill 

• \\Y\\LP(B 1+e \Bi) < 1^ IUp(Bi +e \Bl) + Ce P||y||LP(aBi)- 

Proof of Proposition 3.2. Consider the total decomposition \i of the current 
Ty associated to V. In order to prove Proposition 3.2, we proceed in two 
steps. In the first one (see Section 2.3.1), we apply some elementary operations 
Y |bi +e \Bi j obtaining a new vector field Y\ such that 

• Y\ -= Xb 1+c \b 1 Yi + (xb r - Xsi +e \Si)y still belongs to L\{B R ), 

• \\Yi\\ L p(B r ) < \ \Y\\lp{B r ), 

• divFx = in the interior of B x+e \ B\ , 

• hy 1 \-S — /-ty| Bi XBi \-S where the Borel (for the weak topology) set S 
consists of the 1 -currents R having boundary on dB x+e U dB\. 

In the second step, we modify the currents R £ S (that up to now were 
untouched by our construction). We apply an elementary operation in which 
we cancel (i.e. multiply by 0) the i?'s with both boundaries on dBi +t , and we 
let the others unchanged. Then we consider (identifying the current T 2 with 
a vectorfield Y 2 ) 

T 2 = Y 2 := [ Rd^ Yl (R), 
Js' 

where S' are the currents corresponding to lipschitz curves with one end on 
dBi +e and the other one on dB\ . It follows Y 2 di Y\ and we see that Y 2 is an 
LP -vector field, and since /iy 2 totally decomposes Y 2 , there holds divY" 2 = on 
B 1+e \ B\ . The elementary operations decrease the LP -norms of the boundary 
values, thus 

/ \Y 2 \ p dU 2 < [ \Y\ p dU 2 . (3.1) 

J dBi JdBi 
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We now are in a situation where on one hand 

dT 2 = {v ■ Y 2 \ dBl+e )H 2 LdB 1+£ - (u ■ Y 2 \ dBl )n 2 \-dB 1 , 

where v is the radial vector. On the other hand, by the zero flux condition on 
Y, 

dT 2 LdB 1+e (l) = = dT 2 L8B 1 (l), 

and by homological reasons this implies that the two boundary parts above 
are themselves boundaries. So our strategy is to find another LP -vector field 



Figure 6. We represent schematically, to the left the decom- 
position of (the current associated to) the vector field Y near 
Bi +e \ Bi , in the center the similar decomposition for Yi , and 
to the right the vector field Y 2 , where the part of the decompo- 
sition that will stay unmodified (and does not contribute to Y 2 ) 
is dotted. The result of subtracting Y 2 and adding Y 3 to Y\ can 
be rephrased in a more picturesque way by saying that we are 
"cancelling" Y 2 and "replacing it" by Y3. We eventually loose 
a bit in our estimates, since Y3 "forgets about the support" of 
Y 2 , and no easy form of a superposition principle holds for our 
range of exponents p. 

Y3 whose associated current T3 has dT 3 = —dT 2 \_dBi +e , and which has good 
norm estimates. The choice to which we are led is as follows: 

Y 3 = V/, 

for / solving 

( Af = on B 1+e \ B u 

I dJ = gondB h (3.2) 
[ d v f = on &Bi +e , 

for g := —Y 2 ■ v . Then we can define Y by extending Y3 + Y\ — Y 2 as zero 
outside Bi +e and as Y inside B\ . 

The boundary of the associated current Ty is equal to ( dTy ) Lint^i , therefore 
YeTZ^Bn). 

The only fact left to prove in order to obtain Proposition 3.2, is the estimate 
of the LP -energy of Y , for which we need the following scaling lemma: 



Lemma 3.3. There exists a constant C depending only on the exponent p 
but not on e, such that the following holds. For any f e W 1,p (Bi +€ \ B\) 



REGULARITY OF L p -VECTOR FIELDS WITH INTEGER FLUXES 



23 



that is a weak solution of the Neumann boundary value equation (3.2), where 
g £ L p (dBi) , the following estimate holds: 

||V/||iP( Bl+e \Bi) < Ce~v \\g\\ L p(dB 1 )- 

Proof. We denote by f e a solution of (3.2) with parameter e. We observe 
that the weak formulation of the above Neumann problem states that for all 

( f>eC°°(B 1+e \B 1 ), 

V/ e • V0 = / ct>gdU 2 . 

B 1+e \Bj JdB 1 

Therefore, for any e > and for any test function on IR 3 there holds 

/ V/ £ -V0= / V/i • V0. (3.3) 

JB 1+e \B 1 JB 2 \B 1 

Now observe that the gradients form a closed subspace of L p (B 1+e \ Bi,M. 3 ), 
thus the following equality holds, where q = : 



||V/ £ |U P(Bl+ABl) = sup / V/ £ ■ V0 e : ||V0 e |U, (Bl+ABl) < 1 . (3.4) 

Ubi+A^i J 

Here it is enough to consider functions e belonging to C°°(Bi +e \Bi). For 
any such test function, there exists another test function defined via the 
relation 

e (l + r,0) = 0(l + r/e,0), Vr £ [0,e],V9 e S 2 . 

The map >->■ <p t is bijective between C ca (B 2 \ B t ) and C°°(.Bi+ e \ B t ) and 
for a geometric constant C g < 2 , there holds 

||V&||l,(b 1+ .\Bi) < C s e^ _1 ||V0|| L9( B 2 \ Bl ). 

This last fact and (3.3) can be applied to the equivalent definition (3.4), im- 
mediately yielding our thesis. Indeed, we can obtain a constant C as in the 
theorem's formulation, which depends on C g and on the constant of the clas- 
sical L p -regularity estimate for the Neumann problem on the domain B 2 \Bi, 
neither of which depends on e. □ 

— 1 1 1 1 1 

We thus obtained an estimate of | \Y 3 \ |lp(Bi +e \Bi) via e p \ \Y: 2\ \lp{3b 1 ) > an d 
this suffices because of (3.1). Moreover, — ^2||xp(Bi +e \Bi) < ||^IUp(Bi +e \Bi), 
because {Y 1 — Y 2 )\b 1+£ \b 1 di Y\b 1+ c\Bi- This concludes the proof of Proposition 
3.2. □ 

3.2. Proof of Theorem 3.1. In the proof of Theorem 3.1 it is enough to 
consider the case r = 1 , and suppose B = B R , R > 1 , since the general case 
follows via the scaling of the energy. 

If the Xk and the X would be in L p fl7?. 00 (-B) , then we would apply Propo- 
sition 3.2 to Yi,. = Xk — X on the shell B\ +eo \B e . In general we cannot rely 
on this hypothesis, so we use the fact that IZoo(B) is dense in L%(B) and 
complicate a bit our constructions. 
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Proof of Theorem 3.1: We proceed in 3 steps. 

Step 1: finding a spherical shell of small norm. X G L P (B) and the X k 
converge weakly to it, so by weak lower semincontinuity, up to forgetting the 
first terms of the sequence Xk , there holds 



\Xt — X\\lp < \\X k \\lp + \\X \\lp < 3 1 



LP- 



We fix eo and we divide the interval [1,1 + eo] in M smaller intervals Ih of 
length at least e = eo/2M. Then, with the notation 

A Ih = (i 6 I 3 : \x\ G h}, 

by pigeonhole principle we can find a subsequence of the Xk and an index h 
such that 



\Xk — X 



\p < — 



2C 



From now on we forget about h, and call I := Ih- Given any 5 > 0, up to 
choosing another subsequence and changing I slightly, we can also assume 



\Xk — X\\ p r 



LP(dB iniI ] 



< 5. 



Step 2: approximating the interpolant. At this point, with the notation 
Yk := Xk — X , we use the strong density of IZ^B) in L%(B) to find an 
approximant Yk G IZ^B) such that the L p -distance of Yk and of Yk on Aj, 
as well as the //-distance of their boundary values, are not larger than e\. 
Similarly we can define approximants Xk,X. 

Up to changing I slightly, we can insure that none of the Yk have any charges 
on dAi, so that we can apply Proposition 3.2 to them. We obtain Yk G IZoo(B) 
that is 

• L p -close to Yk on B\, 

• zero om B \ Bi +eo . 

Up to passing to a subsequence there holds: 



X k -Y k ^X k eL p z {B), 
(X k -Y k ) X k \ Bl+eo \ Bl , 

(X k -Y k ) -+X\ Bl . 



The X k defined as above (which depends of the choices of subsequences, on J, 
and on the parameters e±, eo, e, 6), will be our choice of an interpolant between 
X and Xk- 



Step 3: final norm estimates. We can now patch together all our con- 
structions and estimates to obtain the following chain of inequalities. We 
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simplify the notations and write directly || ■ \ \x instead of | 
||^c||i?i < ||^fc||Bi +£0 < ||Xjfc|| Sl+eo by minimality of X k 



\Lp(X) 



< | \X\ \b± + liminf ||ifc||si +eo \Bi by lowersemicontinuity 

< + L71iminf (\\Y k \\ Al + e x + Ce" 1 ] |F fc | | 5B . nf / + e x ) using Prop. 3.2 

< mBi + c (£. + l. 

and since there is no obstruction to letting e, 5 be arbitrarily small, the desired 
inequality 

||^A;IUp(Bi) < H^IUp(Bi), 

holds and the thesis follows. □ 

4. The regularity result 
4.1. Dimension of the singular set. 

Definition 4.1. For a vector field X £ L P (Q) defined on some domain Q, 
we define the regular set of X , reg(X) G Q, as the set of those points in a 
neighborhood of which X is C 1 -regular. The set Q \ reg(X) := sing(X) is 
called the singular set of X . 

Proposition 4.2. // X £ is a minimizer of the L p -energy, then for 

VL' d Q, "H 3_2p (sing(X) R Q') = and sing(X) is nowhere dense in Q' . 

Proof. Without loss of generality we suppose that X is minimizing with respect 
to perturbations supported in a neighborhood N of Q , and we prove the result 
with Q instead of Q' . From Proposition 2.1 we know that xq £ reg(X) if for 
some r > there holds 

r 2p ~ 3 [ \X\ P < e . 

JB{xo,r) 

We can then cover sine; (X) by 25-balls B 25 ,...,B? S contained in N such 
that the balls B^ , having the same centers and radius 5, are disjoint. Now, 
by monotonicity we obtain 



e < S zp ^ / \X\ P , k = l,...,l. 
Jb 

and summing this on k we obtain 



1 f .„._ II* 



l 5 z-*p < _ / \ X \ P < " llLP(n) . (4.1) 

^0 JuB* e 

After choosing such a family of balls for all 5 we obtain the volume estimate 

V? (\Jb{\ = 15 3 < C5 2 p^0, 
therefore by dominated convergence, 

\X\ P -> as 5 ->■ 0. 
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Inserting this in (4.1) gives, by definition of T-L 3 ~ 2p and by the covering property 
of our chosen balls, % 3_2p (sing(X)) = 0, as wanted. 

If we choose a ball B C Vt and we pack it as above with families of small 
disjoint balls of radiuses 5 — > 0, we see by the scaling reasoning as above that 
if X has rescaled energy bounded from below by eo on all balls for all 6 , then 
X has to have infinite energy on B , which is not the case. Therefore there is a 
small ball on which the e-regularity theorem 2.1 holds, showing that sing(X) 
is nowhere dense. □ 

4.2. Singular set of weak limits of minimzing vector fields. 

Proposition 4.3. Suppose that X^ are minimizers and X^ — 1 X . Then 
Xk — > X locally uniformly on Q' \ So, for any Q' (£ Q. Moreover So is 
contained in the energy concentration set 

£ := \ x E Q : liminf limr 2p ~ 3 / \X k \ p > 4e l , 

I k ^°° JB{x,r) J 

where eo is the constant of the e-regularity theorem 2.1, and 'H 3 ^ 2p (T,r\Q') = 0. 

Remark 4.4. It can be shown that S$ = S , but we don't need this character- 
ization. 

Proof. We can assume up to taking a subsequence that X^ — > X strongly 
in L p . We show that T-L^ 2p {Jl) = 0, and that outside £ the X k converge 
uniformly; this is equivalent to the thesis. 

It follows, directly from its definition, that £ can be covered by finitely 
many balls Bi, with centers in £ and radiuses r«, and such that for k large 
enough, 



(2n) 2p - 3 f |X fc | p >2e , for all k,i. 

J2Bi 



>2Bi 

We fix the choice of this set of balls, such that 

i 

Then, by the estimates of the e-regularity, it follows that X k are uniformly 
Holder on Q' \ [J Bi , and therefore they have a subsequence converging uni- 
formly on that set. By the reasoning of the proof of Proposition 4.2, as 5 — > 
the sum J^r 3-2p must converge to zero, and by the arbitrarity of e above it 
follows that H 3 ~ 2p (T) =0. □ 

Corollary 4.5. Let X k be a minimizer of the L p -energy, X k — ^ Xq and 
Sk '■= sing(Xfc) for i > 0, and s > 0. Then for any Q' (s Q there holds 

HlciS n n') > iimsu P ^(5 fc n n'). 

Proof. Consider the balls B k as in Proposition 4.3, except that this time they 
are used to approximate T-L^S) . Then for k large enough there holds 



S k c\J Bi 
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and therefore we can obtain 

n^s nn') + e> lim u^s* n n'), 

fc— >oo 

as wanted. □ 

4.3. Monotonicity and tangent maps. We consider now a sequence of 
blow-ups of a minimzer X around a point Xq. We call X r (x) = ^X{rx + xq) , 
and we observe that 

l6L p z (B r (i„))^I r 6^(B) 

Proposition 4.6. [Monotonicity formula] If X £ L| is a minimizer of the 
L p -energy, then for all x £ B and for almost all r < dist(x, dB) there holds 



d 
dr 



(r 2p ~ 3 [ \X\ p dH 3 ) = 2pr 2p ~ 3 [ \X\ p ~ 2 \X^\ 2 dH 2 (4.2) 

V JB r {x) J JdB r (x) 

where X" is the component of X orthogonal to d/dr. 

Since the right hand side is positive, the left hand side has a limit L(x) for 
r — > + , so we can integrate equation (4.2) from to A, getting 

A 2 P -a I \x\ p dU 3 - L = 2p ( r 2p - 3 \X\ p - 2 \X^ 2 dH 3 . 
Jb x Jb x 

As in [HrL], the function L(x) is actually upper semi-continuous. 
The equation (4.2) also implies that 

/ \X r \ p = E r (X) := r 2p ~ 3 [ \X\ P 

JBi J B r 

is increasing in r, therefore the X\ have a LP -weakly convergent subsequence 
X x . — Xq £ L p , Aj — > 0. By a change of variables in the integrated formula 
we obtain 

A 2p ~ 3 I \X\ p dU 3 - L = 2p ( r 2p - 3 |X A | p - 2 |xj[| 2 ^ 3 , 

J B x JB X 

therefore 

lim / r 2p - 3 \X x \ p - 2 \xl\ 2 dH 3 = 0, (4.3) 

A ~> 0+ J Bi 

Since p' = -£j and X £ L p , we obtain that |X | P ~ 2 X £ // ; the weight r 2p ~ 3 
actually worsens the convergence above since it's bounded away from zero, so 
we obtain that Xq = . This proves more in general the following: 

Proposition 4.7. For any minimizer X , for any x £ int(S) and for any 

sequence of rescalings X x \. around x, with \ — > 0, the weak accumulation 
points X X)Q are radially directed. 
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4.4. Stationarity and dimension reduction for the singular set. From 
now on we call s any exponent (smaller than 3 — 2p, as seen above) for which 
"H s (iS n Q') > 0, where S = sing(X) for a minimizer X . Except for xq in a 
set 5" such that "H S (S") = 0, there holds 

lim inf X~ S K S (S n B x/2 ) > 0, (4.4) 

where the balls Bx/2 are all centered at x . As in the previous section, for 
a subsequence Aj — » our blow-ups converge to a radial tangent map Xq 
weakly in LP , and since they are all minimizers, by Theorem 3.1 they converge 
strongly, up to taking another subsequence, and also Xq is a minimizer. 

The singular set Si of X Xi is the blowup of S , and 
\r s H s (S D B Xi/2 ) =H s (Sin B 1/2 ) 
and from (4.4) we follow that 

H a (S Q nB 1/2 )>0, (4.5) 
where So is the singular set of X . 

Using the radial direction of Xq and the stationarity (Prop. 5.1) we now 
show the following fact. 

Lemma 4.8. For any minimizer X , any tangent map Xq satisfies 

\Xq\{x) = \x\- 2 \Xq\{x/\x\). 

Proof. We use the equation (5.1) with respect to a local frame ei,e 2 ,e 3 such 
that the vector e 3 is the radial one and u associated to X has just the compo- 
nent parallel to de 1 Ade 2 different from zero (as was proved in Proposition 4.7), 
and we consider a perturbation field that can be expressed in polar coordinates 
(p,9) as 

v(p,o) = f(p)4>(o)p. 

We then get from (5.1) that 

= pj \oon P ,e)- p f(p)p 2 dp ( f ) (e)d0- J \u\v(p 1 6)j 2 d p (p 2 fXp))p 2 dpci ) (6)de, 

By the arbitrarily of <f>(9) this translates into the following equation holding 
for almost all 9 

f M p (p, 9) [2(p - l)pf(p) - p 2 f'(p)] dp = 0. 
This can also be written in terms of F(p) = p~ 2p f(p)] as 

\unp,9)p 2 ^F'(p)dp = 



and since this holds for all F with support contained in ]0,oo[, it must be 
that 

\co\(p, 9)p 2 is independent of p, 
as wanted. □ 
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Along the same lines as the above proof (we just have to redefine the the 
orthonormal frames properly), we obtain the following result without difficulty: 

Lemma 4.9. If X\ is parallel to one coordinate direction e% and if the sta- 
tionarity equation holds, then Xi is almost everywhere independent of the 
coordinate x$ . In particular the thesis is satisfied if X\ minimizes the energy. 

Remark 4.10. We note that in Sections 4-3 and 4-4 until this point just the 
monotonicity and stationarity formulas were used, without the intervention of 
any comparison argument. Thus the results proved so far in this subsection 
are valid not only when X is a minimizer, but also when X is just stationary, 
i.e. the 2 -form F associated to it satisfies 



for all families of diffeomorphisms (j) t : B 3 — > B 3 that are differentiably depen- 
dent on t G [—1, 1], equal to the identity in a neighborhood of dB 3 , and such 
that 0o = idffA ■ This requirement is indeed enough to prove stationarity and 
monotonicity. On the contrary, the dimension reduction technique that we are 
about to prove uses the strong convergence result which in turn depends on a 
comparison argument, thus the following proofs hold only for minimizers X . 
An intriguing open question is whether or not the uniqueness of tangent maps 
holds in our case (see Section 6.2 of [Pe3] for a broader discussion) . 

We are now ready to apply the dimension reduction technique of Federer 
to our minimizing vector field X. We start with a radial tangent map X , 
obtained by blow-up at a point x at which So has positive density with 
respect to H s for some s < 3 — 2p as above, as in (4.5). 
As we saw in Section 4.3, Xq is a strong limit of a blowup sequence relative 
to some Xi — )■ + . We also know that the singular set So of X has zero 
"H 3 ~ 2p -measure and is nowhere dense. It follows from Lemma 4.8, that |Xo| 
must be (—2) -homogeneous, and divXo = locally outside Sq. Therefore Xq 
is itself (—2) -homogeneous outside So, and Sq is radially invariant, i.e. 



Now we show that So = {0}. Indeed, were this not the case, we could find a 
point x\ G So fl B 1 / 2 . In this case we could blow up again X with center xi, 
obtaining a tangent map X\ . By strong convergence we obtain that X 1 would 
have to be both directed radially and directed along one fixed direction: this 
would imply that X\ — , contradicting the fact that X\ G Sq ■ 

The following proposition summarizes the above discussion. 

Proposition 4.11. For a minimizing vector field X , the singular set of any 
tangent map sing(Xo) is either empty or contains just the origin. 

After Proposition 4.11 we deduce our main result easily. 

Theorem 4.12. A minimizer X must have finitely many isolated singularities 
in any open Q! d Q. 




XSq d Si 



VA > 0. 
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Proof. If X had an accumulating sequence of singular points sing(X) 3 Xi — > 
x G Q' , then we can select a small r > such that B(x, r) C Q' . Then we can 
consider the distances 

\x X{ | 

K = — ~ A — , 

and we observe that for the blowups Si of ratio Aj and center x, there holds 
H°(SinB 1/2 ) > 2. This contradicts Proposition 4.11 (where H°(S ) < 1) and 
the semicontinuity proved in Corollary 4.5. □ 



5. Stationarity and monotonicity 

5.1. Stationarity formula. We consider a smooth diffeomorphism tp t := id+ 
tV , where V is a compactly supported vector field and t is small enough. We 
compute the stationarity formula arising from 



d 
dt 



\<pM 



= 0. 

t=o 



We recall the formula of the norm of the pullback of uj via tp t , with respect to 
an orthonormal frame field ei, e$: 



\(<p* t u)) x \ 2 = ^^^{dipteiidipte^' 



*ij'=i 

n 

Ei 1 2 

|^VO)( e « + * dV ■ e u ej + t dV ■ e^-) | . 

To deal better with the above t -derivative, we change variable (we let y := 
ift l (x)), so that the point at which we calculate the norm of u does not depend 
on t: 



J \(f* t u x \ p dx = J 



p/2 

^2 K( e * +tdV ■ ei, ej +tdV ■ e j ) | 2 ) det(id + tdV)" 1 



Now we take the derivative of the integrand in t = 0, obtaining by easy 
computations (see for example [Pr]): 



p I \u\ p ~ 2 u(e h e 3 )uj{V e y, e,j) - f \uj\ p diW = 0. (5.1) 
J i,j=i J 

The above formula is justified for minimization problems in LP , because we 
are sure that the manipulations done extend to that setting. What ensures 
that doing the pullback preserves the property of being in L% as well, is the 
following: 

Proposition 5.1. Consider a regular foliation 

{El: A €[-6,6]}, 



REGULARITY OF L p -VECTOR FIELDS WITH INTEGER FLUXES 31 

i.e. a parameterized set of 2 -surfaces in K 3 such that if N e T, = then 
the following (has sense and) holds: 

X-v^dU 3 - [ [ X-u^dH 2 d\ 

where z/ E 2 is the normal vector of S| . 

The following property is equivalent to the fact that X G •' 
For almost all X G [— e, e] the following holds: 



X ■ v dU 2 G Z. (5.2) 

Proof. This follows since TZ°° is dense in L\ in the L p -norm. Suppose indeed 
that there exists a closed C 2 -surface S such that for a set of A G [— e, e] of 
measure 5 > there holds 

/ X ■ vdH 2 E]a + c, a + 1 — c[, for some a G Z. 

In particular, whenever Xi — > X , Xi G then 

\Xi-X\ p dH 3 > C f [ \Xi - X\ v dU 2 dX 



J 



> c I |SaI 1_p 



X.-vdU 2 - / X-vdU 



V 

dX 



> CSc p , 

contradicting the convergence in L p -norm stated above. □ 

Since for t < \\X\\oo/2, it follows that <pt is a diffeomorphism and the 
integral of u on a sphere S is by definition the same as the integral of <f>* t i>j 
on i " 1 (S'), we see by the above proposition that to G ^f(^) implies that also 
the perturbations cj)* t uj belong to the same space for t small. 

5.2. Monotonicity formula. In this section we prove a refinement of the 
stationarity formula. Since the proof is independent if the dimension n of our 
domain, we give a formulation in any dimension (the defnition of J^ifl) now 
requiring the degree to be an integer on any 2 -dimensional sphere). For our 
applications we will just use the case n = 3. 

Proposition 5.2 (Monotonicity formula). If u £ is stationary, then for 
all x and almost all r G]0,i?] with the constraint Br(x) C Vt there holds 



d 
dr 



r 2p ~ n I \u\ p dyj = 2pr 2p ~ n / \u\ p - 2 \d p ju\ 2 da (5.3) 



where d p = ^- is the radial derivative. 
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Proof. We use a strategy similar to [HrL] . If F : Br — > is a weakly differen- 
tiable bijective Lipschitz function, and if u £ then also £ so it is a 

competitor in our minimization. Therefore the stationarity ^ |F t *o;| p = 

0, holds provided that F = ids R and that the family F t is differentiate in 
t. Such properties will be clear from our choices of the map F . (5.3) follows 
from this. 

Definition of F 

Fix < r < s < R such that < t < s/r. Then we define a function 
F = F r s>t : B R — > B R by F(x) := r](\x\)x, such that 

p := \x\ I—)- |-F(x)| 

is continuous and affine on each of the intervals [0, r], [r, s], [s, i?] . We define 

i ;[pi[ S , r] (5-4) 

and 77|] riS [ is defined accordingly: 

s — tr 1 rs(t — 1) 

^Im(p) := - — - + 



s — r p s — r 

Expression of |F*u;| 2 
We do our computation in coordinates. We choose a basis {eo, e±, . . . , e n _i} 
with respect to which to write the matrix dF x , where eo = d p and the other 
vectors form an orthogonal basis together with it. Then 

dF 

— = 7]e k + pr]'5 0k e . (5.5) 

Then 

\{F*u) x \ 2 = [u F (x){dF x ej,dF x e 

(dF dF 



i.j 

n-l 



\ dxi ' dxj 



£ 

tj=0 

i,j>0 i>0 

v^^l + ^iv + PvTld, ■■•' 2 

i,j>0 



The derivative in t 

We now start the computations for the monotonicity formula. 

\F*u\ p = I + II + III (5.6) 

R 
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where, after a change of variables y = F~ 1 (x) , 



I : = / \F*u\ p = t 2p ~ n I \cu\ p dy, 



II : 
III : 



\F*co 



B s \B r 



B R \B S 



*, ,\P 



\F*u 



\uj\ p dy 



Br\B s 



We want now to change variable also in II and to take 4 L , of the terms 

° at I t=l 

above. The easy terms give: 



j/ d 



t=1 (I) = (2p-n) I \oo\ p dy + r / \u\ p da 

B r J QB r 



Ingredients for the computations 



We observe that 



v(p) + pv'(p) 



s — tr 



s — r 



which has t-derivative . It is useful to keep in mind that 77 = 1 for 
t — 1 ; this will be used without mention in the calculations. 
If y = F(x) and a := \y\, then we can write the expression of rj in 
terms of a: 

s — tv 

a = prj(p) = (p — r) h tr 



s — r 



so 



and 



P = f( a ) ■= ( s ~ r )~ — — 4 



v(M) = — + 



s — r 



s — tr 



a -tr 

(s — r) h r 

s — tr 



-1 rs{t- 1) 



s — r 



whence 



r rs 
+ 



s — r a(s — r) 
From (5.5) it follows that for p := |x| G [r, s], 

J(dF- l )=[r ! ( P r- 1 (r ! (p)+ P r ! '(p))Y 1 



so 



dt 



J(dF~ l ) 



t=\ 



(1 ~ n) dt" 1 



+ 



d ( s — r 



:i-n) 



t=l 
r 



dt \s — tr 
rs 



t=i 



s — r <t(s — r) 



+ 



s — r 
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The computation of the i-derivative 

We call I^Wy] 2 = J2ij>o u ij an d we obtain 



II 
II' 

(derivative of the domain) 
(der. of the Jacobian) 
(der. of the main term) 



= r 


1 


|i*w y |V + 2 


J B s \B r t 


d 

dt 


t=i 


= —r 




\ujy \ p da 


J< 





s — tr 



s — r 



p/2 



V 2 \u y (y,-)\ 2 ) J{dF- l )dy 



+(n-l) 



V 
2 



s — r 



B s \B r 



1 - — \ou y \ p dy + 

\y\ I s — r 



\u! y \ p dy 



B s \B r 



A\u\ p [ — -l)dy 



s — r 



_s -r J Bs \B r 

We now take the limit s \, r and we are interested in seeing what the 
equation I' + IV + III' = becomes. The answer is 



/ A\u\ p - 2 \d p ^co\ 2 dy ■ 

J B s \B r 



lim //' 



-r 

IdB-, 

+0 + r 



u\ p da 
\u\ p dy 



dB r 



+0-2pr / \u\ p ~ 2 \d p ju\ 2 dy 

J dB r 

-2pr [ \uj\ p - 2 \d p ju\ 2 dy, 

JdB r 



and 



lim/' = (2p — n) I \co\ p dy + r I 

S ^ T J B r JdBr 



\u\ p dy 



Summing up and using the fact that a; is a minimizer of the energy, we get 
(2p-n) / \u\ p dy + r / \u\ p dy = 2pr / \u\ p - 2 \d p ju\ 2 dy 

J B r J QB r J dB r 

Multiplying both the r.h.s. and the l.h.s of the above equation by r 
get the desired formula 

4- ( r 2p ~ n f \u\ p dy] = 2pr 2p - n [ \u\ p - 2 \d p ^u\ 2 dy 

dr V JB r J JdBr 



□ 

In terms of vector fields, we can state the following: 

Proposition 5.3 (Monotonicity formula, alternative formulation). // X G L% 

minimzes the energy, then for almost all r e [0, R] there holds 

d 



dr 



f r 2 P -n f \ X \ p dy] =2pr 2p ~ n [ \X\ P ~ 2 \X - (X,v Br )v Br \ 2 dU 2 (5.7) 

V Jb t J JdBr 
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